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Abstract. We investigate the effect of hole and electron doping to half- filling in the periodic 
^ , Anderson model on a triangular lattice by the Hartree-Fock approximation at zero temperature. 

$_^ ' At half-filling, the system exhibits a partially disordered insulating state, in which a coUinear 

jy^ ' antiferromagnetic order on an unfrustrated honeycomb subnetwork coexists with nonmagnetic 

state at the remaining sites. We find that the carrier doping destabilizes the partially disordered 
C^ ' state, resulting in a phase separation to a doped metallic state with different magnetic order. 

The partially disordered state is restricted to the half-filled insulating case, while its metallic 

counterpart is obtained as a metastable state in a narrow electron doped region. 
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1. Introduction 

Partial disorder is one of the peculiar orderings to geometrically-frustrated systems. It is 
stabilized by relieving the frustration with self-organizing the system into coexistence of 
^^ . magnetically-ordered sites and nonmagnetic sites. Interestingly, the partially disordered (PD) 

^T I state has been observed in some metallic compounds [1, 2j, for example, in Laves compounds 

RMn2 (i?=Dy, Tb, Th) and in /-electron compounds such as UNi4B [3,4] and CePdAl [5j_6j. In 
the Laves compounds, the PD state appears for Mn spins as a mixed state of ordered moments 
and nonmagnetic sites in the transient regime between magnetically-ordered and paramagnetic 
states while changing the lattice parameter. On the other hand, in the /-electron compounds, 
the PD state shows a character of a mixture of ordered moments and Kondo spin singlets. For 
this phenomenon, therefore, an instability toward Kondo singlet formation plays an important 
role. While the coexistence of ordered moments and spin singlets is expected to cause interesting 
C^ ' magnetotransport phenomena, the comprehensive understanding of the PD state is not reached 

yet both experimentally and theoretically. 

In the present study, we focus on such a PD state with coexistence of magnetic ordering 
and nonmagnetic sites related with Kondo singlet formation. Such coexistence was theoretically 
studied as a consequence of competition between the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction O El [9] and the Kondo coupling [lOl Hi] under geometrical frustration. The PD 
state was indeed obtained by the mean-field (MF) approximation for a classical pseudomoment 
model jf2| [T3] and by a sophisticated variational Monte Carlo method for the Kondo lattice 
systems [13] • Recently, the authors found a similar PD state in the periodic Anderson model 
on a triangular lattice at half-filling by the Hartree-Fock approximation [15j. The PD state, 
however, is insulating, while it is metallic experimentally. 
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In this paper, to search for a metalhc PD state, we extend the analysis of the PD state for 
the periodic Anderson model by hole and electron doping to half-filling. We will demonstrate 
that the PD state is immediately destabilized by carrier doping at the MF level, with showing 
a phase separation to a doped magnetically-ordered state. We show that a metallic PD state is 
obtained only as a metastable state in the MF solutions for electron doping. 

2. Model and method 

We consider a periodic Anderson model on a triangular lattice, whose Hamiltonian is given by 

^ = -* E (4.C,. + H.c.) - V 5^(4/,. + H.c.) + [/ ^ 4/.t/i/4 + E,Y. /-/-' (1) 

(i,j),a- i,cr i i,a 

where c]^ (Qo-) and //^ {fia) are the creation (annihilation) operators of conduction and localized 
electrons at site i and spin a, respectively. The sum of {i,j) is taken over the nearest-neighbor 
sites on the triangular lattice. Hereafter, we take the hopping t = 1 as an energy unit and the 
c-f hybridization V > 0. In the following calculations, we fix the Coulomb interaction between 
localized electrons U = 2 and the atomic energy of localized electrons Eq = —1. 

In the previous study [15], the authors found a PD state at half-filling in the ground state 
of the model ([1]) by applying the Hartree-Fock approximation to decouple the U term with 
considering three-site unit cell. The PD state is stabilized between a 120° antiferromagnetic 
(AF) metal and a Kondo insulator (KI). Here, to explore a metallic PD state, we extend the 
analysis to the doped regions near half-filling by employing the same method as in Ref. [15]. For 
better resolution, we take the sum over 400 x 400 grid points in the three-folded Brillouin zone 
in calculating the mean fields in this paper. 

3. Result and Discussion 

Figure 1 shows the result of the ground-state phase diagram including the carrier-doped regions 
near half- filling at U = 2 and Eq = —1. Schematic figures for spin and charge states in each 
phase are also shown, n is the electron density defined by n = ^^ ^(q ^q^o- + fiafi,^) /-^ ^ where 
N is the total number of sites. At half-filling n = 2, a sequence of three phases — an AF metal, 
PD insulator, and KI — is obtained while increasing F, as in the previous study [15j . In the 
PD state in the intermediate V region, one of three sublattices becomes nonmagnetic, while a 
collinear AF order appears in the remaining two sublattices. 

When we dope carriers, the PD state is unstable to a phase separation (PS) for both electron 
and hole doping, as shown in figure 1. For electron doping, PS takes place between the PD state 
and the ferrimagnetic (F) metallic state. (In the very narrow window of 0.875 < re < 0.883, PS 
takes place to the AF metal.) Meanwhile, for hole doping, PS occurs between the PD state and 
the AF metallic state extended from half-filling in the small V region. 

The regions of PS are determined by the following procedure on the basis of grand canonical 
ensemble. Figure WL.^) shows the zero-temperature grand-canonical potential ^ = E — fin for 
various ordered states measured from that for the F state {E is the internal energy per site and 
fi is the chemical potential). For each /i, the lowest Cl state gives the ground state. Figure [2]^b) 
shows the electron density re as a function of fx. In general, a transition between different ordered 
states is of first order, and hence, re shows a discontinuity. The PS region is identified as the 
region where re jumps. In the present example at V = 0.90, the ground state changes from the 
AF state to the PD state at ;U ~ 0.244, and finally to the F state at /i ~ 0.391, as shown in 
figure[2]^a). Accordingly, re changes discontinuously from 1.93 to 2.00 at /i ~ 0.244 and from 2.00 
to 2.26 at ^ ~ 0.391, as shown in figure [2]J^b). These two regions are PS: the former is between 
the AF metal and PD insulator, while the latter between the PD insulator and F metal. 
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Figure 1. Ground-state phase diagram of 
the model in eq. (1) at C/ = 2 and Eq = 
—1. Schematic pictures for each phase are also 
shown: Upper (lower) circles show the c (/) 
electron. The sizes of the circles represent local 
electron densities, and the arrows represent local 
spin moments, schematically. A hatched area 
shows the parameter region where a metallic PD 
state is obtained as a metastable state. 



The PS regions are also identified by the canonical ensemble method. Figure [5][|c) shows the 
internal energy i? as a function of n. In this approach, the state with the lowest E gives the 
ground state for each n, and the PS regions are determined by drawing tangent lines for E 
curves for different ordered states. Kt V = 0.90, as shown in figure [2I|c) , two tangent lines can 
be drawn from the PD state at half-filling n = 2: one is to the AF state at n ~ 1.93 and the 
other is to the F state at n ~ 2.26. These indicate two PS regions, consistent with the above 
results by the grand canonical ensemble scheme. 

We conducted similar analyses while varying V to complete the phase diagram shown in 
figured! As a result, we found that the PD state is destabilized by carrier doping to half- filling 
in the parameter range considered here. That is to say, the PD state is limited to the insulating 
state at half-filling and a metallic PD state does not appear as a ground state within the present 
MF calculation. 

In contrast, the carrier doping to the AF state retains the AF order, and does not lead to PS. 
In particular, for hole doping, the AF metallic state extends widely, as shown in figure [TJ For 
electron doping, however, the AF state terminates at n ~ 2.05, and shows PS to the F state. On 
the other hand, the KI state changes into the F state without showing PS for electron doping, 
whereas it immediately exhibits PS to the AF state for hole doping. 

Although the PD state is unstable against the doping, we note that a doped metallic PD state 
is obtained as a metastable MF solution. The metastable metallic PD solution exists in the small 
region of 0.88 ^V ^ 0.96 for electron doping as indicated by the hatched area in figure [TJ The 
corresponding Vt and E are shown in figure [2)^a) and figure EJc) , respectively. The proximity 
to a metallic PD state can be understood from the electronic structure at half-filling shown in 
figure 3 and figure 4 in Ref. [15]. For 0.88 ^V < 0.96, the lowest unoccupied band (LUB) has 
an energy minimum around the M' point in the Brillouin zone, and has a considerable curvature 
in the band bottom. An electron doping to this band bottom retains the PD state in the doped 
region. On the other hand, in the larger V region, the bottom of the LUB becomes flat with less 
curvature, and accordingly, the density of states (DOS) has a sharp peak at the upper edge of 
the gap. Doping to the flat band leads to an instability toward PS. Similar behaviors are seen 
in the hole doped case: For hole doping, DOS shows a sharp peak at the lower edge of the gap 
in the entire range of V , and hence, a metastable PD metal is not found there. 

Let us finally discuss a possibility to stabilize a metallic PD state on the basis of our results. 
In our results, the PD state is unstable against PS with the AF metal or the F metal. One 
may extend the PD state to doped regions by destabilizing the AF and F states with additional 
perturbation. A possibility is the AF Ising interaction between localized electrons, Iz Ylu j) SfSj, 

where Sf = {fl^fi'\ — fJ\fii)/'2'. This term was shown to stabilize the PD state at half-filling 
in the previous study for the Kondo lattice systems [1^. Another possibility is the Coulomb 
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Figure 2. (a) Grand canonical 
potential at T = for each 
state measured from that for the F 
state as a function of the chemical 
potential fi. (b) Electron density 
n as a function of //. The shaded 
regions show the PS regions, (c) 
Internal energy per site E plotted 
as a function of n. The straight 
black lines show tangent lines to 
determine the PS regions. All the 
data are at U = 2, Eq = —1, and 
V = 0.90. 



interaction between conduction electrons. In particular, a long-range Coulomb interaction may 
prevent the system from PS. It will be interesting to examine whether a PD metal appears when 
PS is suppressed. It is also important to consider other ordered states in larger unit cells beyond 
three-sublattice orderings in the present study. 

4. Summary 

We have investigated the carrier doping to a partially disordered state in the periodic Anderson 
model on a triangular lattice by the Hartree-Fock approximation. Elucidating the ground state 
phase diagram near half-filling, we found that hole and electron doping to the partially disordered 
phase leads to a phase separation to a doped magnetically-ordered state. In the small c-f 
hybridization region, a metallic partially disordered state is obtained as a metastable solution 
for electron doping. We have discussed possibilities to stabilize the metallic partially disordered 
state on the basis of the results. 
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